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Recent measurement of time dependent asymmetries in B — > <pK are indicative of a 
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new source of CP violation. We examine squark mixing in SUSY as this new source, 
Q_i| and using QCD improved factorization method to describe B — > <f>K decay, find the 

Oh. allowed range of parameter space for p and ip, the magnitude and phase of the down 

^ . type LR(RL) squark mixing parameter & lr(rL) ■ then study B — ► (j)K* and calculate 

the expected CP asymmetries in the same range of parameter space. We find that this 
asymmetry is in the range 15% to 20% for acceptable value of B — > <pK* branching 

h ; 

ratio. We also predict the helicity dependent CP asymmetries in the same parameter 
space. 
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1 Introduction 



Time dependent asymmetries measured in the decay B — > <pKs both by BaBar and Belle 
collaborations [1, 2, 3, 4] show significant deviation from the standard model and this has 
generated much theoretical speculation regarding physics beyond the standard model [5-20] . 
In the standard model, the process B — > 4>Ks is purely penguin dominated and the leading 
contribution has no weak phase. The coefficient of sm(Am,Bt) in the asymmetry therefore 
should measure sin 2/?, the same quantity that is involved in B — > ipKs in the standard 
model. The most recent measured average values of asymmetries are [4, 21] 

S i , Ks = 0.734 ±0.055 

S^ Ks = -0.15 ±0.33 (1) 

The value for S^k s agrees with theoretical expectation from the CKM matrix of S^k s = 
sin 2(3 = 0.7151^045 [22]. This leads to the conclusion that CP phase in B — B mixing is 
consistent with the standard model. The deviation in the <pKs is intriguing because a penguin 
process being a loop induced process is particularly sensitive to new physics which can 
manifest itself in a loop diagram through exchange of heavy particles. In this article we will 
consider effects arising from non universal squark mixing in the second and third generation of 
the down type squarks in supersymmetric theory as the origin of additional contributions to 
the amplitude within the mass insertion approximation scheme. In particular, the exchanges 
of gluinos (g) and squark (q) with left-right mixing can enhance the Wilson coefficient of 
the gluonic dipole penguin operator 0$ g by a factor of rrig/mb compared with the standard 
model prediction and we take into account its effect on the process B — > (pKs- In our 
analysis we take the B — B mixing phase the same as in the standard model as required by 
ipK s data, and permitted in SUSY by requiring that the first and third generation squark 
mixing to be small. We study B — > <pK in QCD improved factorization scheme (BBNS 
approach ) [23]. This method incorporates elements of naive factorization approach (as its 
leading term ) and perturbative QCD corrections (as sub-leading contributions) and allows 
one to compute systematic radiative corrections to the naive factorization for the hadronic B 
decays. Recently several studies of B — > PV, and specifically B — > <pK have been performed 
within the frame work of QCD improved factorization scheme [24-28]. In our analysis of 
B — > <f)K, we follow [28] which is based on the original paper [23]. In this analysis, we 
only consider the contribution of leading twist meson wave functions, and also neglect the 
weak annihilation contribution which are expected to be small. Both these introduce more 
model dependence in the calculation through the parameterization of an integral, which is 
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otherwise infrared divergent. 

In supersymmetry, assuming masses of squarks (q) and gluinos (g), the new source of 
CP violation can be parameterized by the complex quantity 5 b £ R ^ RL ^ written in the form 
p e »V> \y e identify the region in p — ip plane allowed by the experimental data on B — > (f)K 
time dependent asymmetries S^Ks and C<t>K s and the branching ratio. This allowed region 
is dependent of the QCD scale /z, therefore we illustrate the region for two values of /i = m b 
and m b /2. The same contribution should also be present in other penguin mediated process. 
We study the effect of LR(RL) mass insertion to the B — > <f)K* decay mode which is also 
a pure penguin process using QCD improved factorization method. We then estimate the 
branching ratio BiB — > (j)K*) and the CP asymmetry Acp in the parameter space of S b ^ R ^ RL ^ 
allowed by B — > (f)K data. In this vector vector final state, one can also construct more CP 
violating observables [30]. We compute these observables in the same range of parameter 
space as that allowed by B — > <pK. 



2 CP Asymmetry of B — > (j)K 

The time dependent CP asymmetry of B — > (j)K s is described by : 

* Ks[) " r(F(^^ s )+r(so ( ^^) [Z) 

= -C^Ks cos(Am B t) + S^ K sm(Am B t) (3) 
where S^k and C^k s are given by 
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1+ I ^4>K S 

and A^^g can be expressed in terms of decay amplitudes: 

_ 2i0 M{W ^ <t>K s ) 
X * Ks ~ " 6 MtB^TWi) { ) 

The branching ratio and the direct CP asymmetries of both the charged and neutral 
modes of B — > <fiK have been measured [1, 2, 3, 4, 21, 31]: 

B(B° -> = (8.0 ± 1.3) x 10~ 6 (6) 

B(B+^(f)K + ) = (9.4 ±0.9) x 10~ 6 , (7) 
S (j}Ks = +0.45 ±0.43 ±0.07 (BaBar); (8) 
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scale 


Ci 


c 2 


c 3 


C 4 


c 5 


c 6 


fj, = m b /2 
H = m b 


1.137 
1.081 


-0.295 
-0.190 


0.021 
0.014 


-0.051 
-0.036 


0.010 
0.009 


-0.065 
-0.042 




CV/ «era 


C$/ a em 


C%j a em 


Clo/ «em 




C 8 g 


fi = m b /2 


-0.024 


0.096 


-1.325 


0.331 


-0.364 


-0.169 


H = m b 


-0.011 


0.060 


-1.254 


0.223 


-0.318 


-0.151 



Table 1: Standard model Wilson coefficients in NDR scheme. 



= -0.96 ± 0.50+°;?? (Belle); (9) 
C 4>Ks = -0.19 ±0.30 (10) 
A CP (B + -> (f)K + ) = (3.9 ±8.8 ±1.1)% (11) 



3 The exclusive B —> <pK decay 

In the standard model, the effective Hamiltonian for charmless B — * <pK{(f)K*) decay is given 
by [23] 

G p 
V2 



Keff = 7^ V tbKs 



10 

(12) 



c^om + c 2 (n)o 2 (n) + J2d(v)Oi(ri + c 7l o 7l + c Sg o 8g 

i=3 

where the Wilson coefficients Ci(fx) are obtained from the weak scale down to scale fx by 
running the renormalization group equations. The definitions of the operators can be found 
in Ref.[23]. The Wilson coefficients C, can be computed using different schemes [32]. In this 
paper we will use the NDR scheme. The NLO values of Ci(i = 1 — 10) and LO values of 
CVp C 8g respectively at /i — m b /2 and m b used by us based on Ref. [23] are shown in Table 1. 



4 B —> (f)K in the QCDF Approach 

In the QCD improved factorization scheme, the B — > <pK decay amplitude due to a particular 
operator can be represented in following form : 

< <j)K I O I B >=< <f>K\0\B > fact [l + £r„< + 0(A QCD /m b )] (13) 

where < (f)K \ O | B >f ac t denotes the naive factorization result. The second and third term 
in the bracket represent higher order a s and AQco/n^b correction to the hadronic transition 
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amplitude. Following the scheme and notations presented in Ref.[28, 29], we write down the 
B — > (j)K amplitude in the heavy quark limit. 

M{B + ^(j)K + ) = 

M(B°^^K°) = ^mlUF^ K {ml)V pb V; s [al + al + al 



( a 7 + a 9 ~j~ a lp) _|_ v 

o r a 10a 



where p is summed over u and c. The coefficients af are given by 

C 4 



(14) 
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(15) 



with CV = (N% — 1)/2JV C and 7V C = 3. The quantities V^, H^, P 1 ^ and are hadronic 
parameters that contain all nonperturbative dynamics. 
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rrih Jo Jo 



(1-x) 



(16) 



where, Sj = mf/m 2 . Here, represent contributions from the vertex correction and 
correspond to hard gluon-exchange interactions with spectator quarks. P% and represent 
QCD penguin contributions. We neglect order a em EW penguin corrections to a^. fB,fx 
are the B and K meson decay constants and Ff^ K denotes the form factor for B — > K 
transitions. & B (z), $ K (x), and 3>^(y) are the 5, K, and meson wave functions respectively. 
In this analysis we take following forms for them [28] 



m B x 
~2u£ 



(17) 



$ B (x) = N B x 2 (1 - x) 2 exp 
®K,<t>( X ) = 6 x (1 - x) 

where, N B is a normalization factor satisfying Jq dx $b(x) = 1, and = 0.4. 

For the sake of completeness, we give the branching ratio for B — > <pK decay channel in 
the rest frame of the B meson. 



BK(B -> 0X) 



8n m 2 B 



M(B^</>K) 



where, r B represents the B meson lifetime and the kinematical factor | P a 

| Pcm |= T^-yVs - (^A- + «^) 2 ] t m l - ( m X - m <f>) 2 } 



(18) 
is written as 
(19) 



5 SUSY gluino contributions to B — > <f>K 

In order to study the new physics contribution to the CP violating phase of amplitude 
M(B — > 0X), we compute the effect of flavor changing contribution to B — > arising 
from q — q — g interactions in supersymmetric theory under the mass insertion approximation 
scheme [33, 34]. In this approximation, the flavor changing contribution is parameterized in 
terms of 5j[ B = A^ B /m 2 , where, A represents the off-diagonal entries of the squark mass 
matrices, m is an average squark mass, A, B = L,R and i,j are the generation indices. 
The LR(RL) mass insertion can enhance the Wilson coefficients C 77 and C &g by a factor 
of rrig/miy compared to the standard model contribution. This leads to a strong limit of 
order 0(1(T 2 ) on the LR(RL) insertions | S b L s R{RL) | from the B(B -> X s ^) [34, 35] while 
the limit on the LL and RR ones is rather mild [34, 35]. Thus, although larger values for 
LL and RR mixings are allowed, when one considers B — > (f>K, the effect of their mixings 
are only significant in the parameter space where the squark and gluino masses are at the 
edge of their experimental constraints [17]. Motivated by this fact, we only concentrate on 
LR(RL) down type squark mixing in hereafter. Thus, the new physics effect is very sensitive 

to 3 b LR(RL)- 
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In general, these contributions LR(RL) can generate gluonic dipole interactions with 
the same as well as opposite chiral structure as the standard model. In our analysis we 
will consider each of them separately. Furthermore we will only consider the gluonic dipole 
moment operator, which is the dominant operator for this process. 

The effective Wilson coefficient for C^ SY obtained in the mass insertion approximation 
is given by for the same chiral structure as the standard model [36, 37] 



with 



G F (V ub V u % + V ch V* s )ml LR ( RL) m b 



G(x) = X [22 - 20a: - 2a: 2 + 16x Info) - x 2 ln(x) + 9 Info)] , (21) 

Oil X J 

where x = m'-fra'- is the ratio of the gluino and squark mass. 

Using the renormalization group equation one can evolve the coefficient C% g USY from the 
high scale rriq to the scale m b relevant for B — > (pK decay [36] 

C!r(m b ) = ?Ar y K~)> ( 22 ) 

with 

77 = (a s (mq)/a s (m t )) 2/21 (a s (m t )/a s (m b ) ) 2/23 (23) 

One can obtain C§ g USY for opposite chirality, by adding one more operator similar to Os g 
with (1 + 75 ) — > (1 — 75) and 5 b £ R — > 5 RL . However, in B — > <pK process, both LR and RL 
contribute with the same sign because B and K parity are both 0~, and the process is parity 
conserving. 

The effective Wilson coefficient C$ f g f is defined as C^ 7 = C Sg + Cjf 5Y '. This effective 
Clg* will contribute to the amplitude M(B — > <f>K) through the function of Equation 
16. C%J depends on the magnitude and phase of the (S^ R ^ RL ^), value of squark mass (m,) 
and the ratio x (= m|/m|). The variation of C%J with x is determined by the function 
I G(x) I as shown in Figure 1. From this Figure, it is clear that SUSY gluino contribution 
to B — > (j)K first increases with increase in x, and then after some value of x — 0.5, it starts 
decreasing asymptotically with further increase in x. 

The different input parameters and their values used in numerical calculation of branching 
ratio and CP asymmetries are summarized in the Appendix. We treat the internal quark 
masses in the loops as constituent masses rather than current quark masses [38] 
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5.1 LR(RL) mixing 

In this section we study the effect of LR(RL) mixing in B — > <pK process. This LR(RL) 
mixing of the down type squark sector can also affect the B — * 7X3 process and B s — B s 
mixing. Hence we need to take into account the limit on LR(RL) mixing parameter S^ R ^ RL ^ 
from the above two experimental data in the present analysis. In the first case, it has been 
shown in Ref. [34] that from the measurement of B{B — > jX s ) one gets | S*£ R , RL \ |< 1.0 x 10~ 2 
and 3.0 x 10~ 2 for x = 0.3 and 4 respectively, with rriq = 500 GeV. It is interesting to note 
that the lower the x value stronger the limit on | 5 b £ R ^ RL ^ | , which can be explained by the x 
dependent behavior of the C^ SY . 

The current experimental data on B s — B s mixing is AM S > 14.4 ps _1 ( at 95% C.L.) 
[39] . We have found that the LR(RL) mixing does not change the value of AM S significantly 
from the standard model prediction in the allowed range of | S^ R ^ RL ^ \. 

In our analysis we consider rriq = 500 GeV and take two values of x = 0.3 and 4.0, which 
will determine the gluino masses. In Figure 2 we show the la allowed region in p — ip plane 
from B — > (pK data on S^k, C<f>K and B. The gray band indicate the parameter space which 
is allowed by S^k- The area outside the two dotted contours is allowed by C^k, while the 
area enclosed by the solid curves is allowed by the B[B — > 4>K S ) measurement. The region 
(marked by Z) in gray band enclosed by the solid curves is the only parameter space left in 
p — ip plane which is allowed by the experimentally measured S^k, C<f>K and B within la. 

The Figures 2(a) and (6) , correspond to contour plots for x = 0.3 and 4.0 respectively 
at the scale p = rrib. For x = 0.3, we get two allowed regions each at positive and negative 
values of the new phase ip. On the other hand for x = 4.0, we get only one allowed region 
which lies at the negative value of tp and at much higher value of p > 2.2 x 10~ 2 . We have 
noticed before that the constraint on LR(RL) mixing parameter from the B(B — > X s '-f) is 
stronger at x = 0.3 compared to the limit at x = 4.0. This behavior is also reflected in the 
B — > <f)K process, where we find that, for x = 4.0, the la constraint from S^k, C^k and B 
is much weaker compared to the constraint shown in Figure 2(a) correspond to x = 0.3. 

Similar allowed regions are shown in Figure 2(c) for a different choice of the QCD scale 
p = rrib/2. One can see that the allowed parameter space does depend on p. In this case, 
both the allowed regions are confined at the positive value of if). For x = 4.0 ( Figure 2(d)), 
there are no allowed regions. From the S^k an d branching ratio contour one can see that 
the allowed region from B — > <pK require some higher value of p which lies beyond B — > X s ^ 
limit. 

Before we conclude this section, we would like to compare our predictions with some of 
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the existing literatures on B — > <f)K process [10, 11, 17]. 

We agree qualitatively with the results of Ref.[ll, 17] in places where we overlap. Similar 
to our approach, both of these analyses were based upon the QCD improved factorization 
scheme. However, there are some quantitative differences between these papers and our 
analysis. For example, we differ in the choice of squark and gluino masses, the authors of 
the above two papers considered degenerate squark and gluino masses, whereas we have 
considered non-degenerate squark-gluino masses. We have fixed the squark mass at 500 
GeV and considered two values of the gluino masses, determined by the parameter x defined 
earlier. Secondly, we have performed our analysis for two values the QCD scale, fx = mb/2 
and m&. Our results depend strongly on the choice of the ratio x and also on the scale 
fx. However, in a broad sense, we do agree that to satisfy B — > <pK data, one requires 

l*&om l~io- 3 -io- 2 . 

In Ref. [10], authors made a detailed investigation of a scenario in which the LR and RR 
operators co-exist. Moreover, because of the large mixing, the calculation was done in the 
mass eigenbasis with more model dependence than ours. It has been shown in this analysis 
that RR insertion (which arises due to a large mixing between s R and b R ) could show sizable 
effect on S^k, but only for very light gluino mass, near the experimental bound. Such a 
large RR mixing also modify AM S significantly which can be observed at the Tevatron Run 
II. In their second case, they have the combination of both large right-right and left-right 
squark mixing (LR + RR). In this case the squark and gluinos could be sufficiently heavy 
to have no significant enhancement of the AM S . 

From our analysis we observe that SUSY leads to a comprehensive understanding of 
B — > <pK s data though in a very limited parameter space of S^ R ^ RL ^ . In rest of the paper we 
now explore the consequence of such LR(RL) mixing of squarks in the B — > 4>K* process. 

6 B —> <f>K* decay 

In this section we will study the effect of LR(RL) mixing of down type squarks to B — > (pK* 
process through the gluonic dipole moment operator C% g . We will study the B — > <pK* 
process by using the QCD improved factorization. Using this method one can compute 
nonfactorizable corrections to the above process in the heavy quark limit. Recently the 
B — > VV process has been computed using QCD improved factorization method [40]. In 
rest of analysis we will follow Ref. [40]. 

The most general Lorentz invariant decay amplitude for the process B — > VV can be 
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Branching ratio 


Data 


Weighted average 


B+ -> (j)K* + 


BaBar ( 12 l + ?o ±15) x 10~ 6 

cleo (io.6±t:9±i:l) x io^ 6 

Belle (9.4 ± 1.1 ±0.7) x 10" 6 


(9.9 ± 1.23) x 10~ 6 


B° -> 


BaBar (ll.l^|±l.l) x lO" 6 
CLEO (ll.5l|;fli;?) x 10~ 6 
Belle (l0±i;|±8:l) x 10~ 6 


(10.6 ± 1.3) x 10~ 6 



Table 2: Experimental data of B -> decays from BaBar [43], CLEO [44] and Belle [45] 
and their weighted average. 



expressed as 

M(B(p B ) -> Vi(ei,pi)V2(e 2 ,P2)) oc e^e^ + bp Bp ,p Bv + ice^p^pi 



(24) 



where the coefficients c correspond to the p-wave amplitude, and a, b to the mixture of s and 
o? wave amplitudes. Using these a, b and c coefficients one can construct the three helicity 
amplitudes: 

1 



H±± 



2m Vl m V2 
a =f m B p cm c 



m B ~ m Vi ~ m V,) a + 2m BPcm b 



(25) 



where p cm is the center of mass momentum of the vector meson in the B rest frame and 
m Vl (m V2 ) is the mass of the vector meson V^V^). These helicity amplitudes H 00 and H±± 
can be related to the spin amplitudes in the transverse basis (Ao,A\\,A±) defined in terms 
of linear polarization of the vector mesons: 



++ 



Ao = Hoo 

- 4 " = T2 (K 

-4x = ^ 



(26) 



The decay rate can be written as 

Pan 



T(B - VM) 



Hqo I ± I H ++ | ± | H. 



(27) 
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Neglecting the annihilation contributions (which are expected to be small) to B — > <f)K*, 
Hqq and H±± are given by: 



#00 = 



G F a n (4>K*)U 
y/2 2m K * 
Gi 



(m% - m 2 K * - mfj (m B + m K *) Af^*(mJ) 



2 „2 



4m B p 



m B + va K > 



K) 



{m B + m K *)A* K \ml)T 2WgPc V 



7«B + 

where, a n (4>K*) = + + ag — (a™ + ag + a" )/2. The effective parameters appearing 
in the helicity amplitudes iJ o an d #±± are given by 

03 - C3+ ^ + 4^iv: C4F 



(28) 



04 = ^4 + -^ + 



C3 CK S C F 



C 3 [F n + G n {s q ) + G»( S6 )] - d (yG n (s u ) + ^G"( Sc ) 



At 



+ (C 4 + C 6 ) £ + l(C 8 + C 10 ) £ ei G n (s t ) + ^C 9 [e,G"( S(? ) - V(s 6 )] + C n effg G r Q 



< 


^ N c 4ttN c K 


12) 


a™ 


C s a 8 C F 












C 9 a s C F 

= Cl0 + T c + ^N- C C * F - 


9^ Ce 



(5 



where \ q = V q bV* q ,,q' — d,s and the superscript n denotes the polarization state of the 
vector mesons, with n = for helicity 00 states and n = ± for helicity ±± states. Note that 
c effg w hich is the combined Wilson coefficient of the SM and SUSY corresponding to the 
gluonic dipole moment operator depends on n, because both odd and even parity transitions 
are allowed. The QCD penguin loop functions G n (s) are give by 

G°(s) = ---ln-^ + 4 f l dx&fAx) f 1 duu(l-u)\n[s-xu(l-u)}, 
3 3 to& Jo 11 Jo 

G ± (s) = ^-^ln-^- + 4 Cdxg^ix) C duuil - u)\n\s - xuil - u)\ (30) 
3 3 nib Jo Jo 

=F - f'dx^^- f 1 duu{l-u)\ -2\n^ + \n[s-xu{l-u)}+ xu{\ u) \ 
2 Jo x Jo V J \ m b 1 v ;j s _ xu{1 _ u) f 

Similarly we also have leading EW penguin-type diagrams induced by the operators 0\ and 
2 : 



C n e = (^G h (s u ) + ^G h (s c 



(31) 
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The dipole operator O g will give a tree-level contribution proportional to 

G ^- 2 f dx m G ^- 2 f dx £M + h lT hf dx £M. (32) 

y Jo x y Jo x 2 2 Jo x 

The vertex correction F n as two contributions, the hard scattering function /" and ffj from 
hard spectator interactions with a hard gluon exchange between the vector meson and the 
spectator quark of the B meson. 

F n = -12 In A - 18 + /" + /" (33) 

with 

/° = £dx$\{x) (^Y^lnx-Zii^ , 

ff = J^dxgt\x){?> l T ^\nx-^y (34) 
The hard spectator interaction is given by (Vi: recoiled meson, V 2 : emitted meson): 

iV c Aln JO Jo 1] Jo f 



f. 



Jo p Jo T] z Jo 



± J B JV2 

" ~ N c ///,;// : 



7V C ///j.// . .A) p Jo { ' C rf 



where 



TUb + 771 Vi 



(m B + m Fl )Ar(0 =F 2m f Cm V BV ^ml 2 ), (36) 

where, the asymptotic form of the leading twist ($^(x), $^(x)) and twist-3 LCDAs (^' 1 '' > 
for the vector meson are defined as [41] 

§V( X ) = $V( X ) = gV> = Q X (1 - x) 



^ 4 



l + (2x-iy , (37) 



and the B meson wave function &b(x) is defined in Section 4. The expression ffj contains 
logarithmic and linear divergent terms arising from the dominance of soft gluon interac- 
tions of the spectator quark. Hence, the factorization breaks down at twist-3 order for 
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X 


IP) V) 


»susy /• 11T1 i+o n f -\r\-G\ 

ID ^111 UI111S Ul 1U J 


•A-CP I 111 /0 J 


0.3 


l U.ft A 1U , -u.o J 


oo 07+4.88 /r,-, 7fi +4.56\ 


— zL 7(— A A\ 


(0.4 x lO" 2 , -0.7) 


21.50ltg (20.17i| 22 ) 


-7.0 (-6.5) 


(0.6 x 10" 2 , 1.5) 


27.46±|;I 5 (26.82±t;g|) 


17.7 (15.7) 


4.0 


(2.4 x 10" 2 ,-0.5) 


24.331^ (22.65±tF) 


_4.7 (-4.4) 


(2.6 x lO" 2 , -0.45) 


26.9811? (25.11+t;^) 


-4.2 (-3.9) 


(2.8 x lO" 2 , -0.8) 


25.31111 (23.8711;^) 


-8.0 (-7.5) 



Table 3: B(B+ -> 0fsT*+) and „4 C p(5 + -> 0^*+) at the QCD scale p = m b for LR mass 
insertion for selected points in the allowed p — tp space. The numbers in the parenthesis 
correspond to the RL mass insertion. The standard model branching ratio corresponding to 
this scale is (6.18±i;?§) x 10 6 . The errors are due to ±10% theoretical uncertainties in the 
calculation. 



transversely polarized vector mesons. We parameterize these divergent integrals in the form 
(1 + pbexpi(f>b) lnrriB/AQC'D, with an arbitrary phase (fib and pb < 1. In our analysis, we take 
Pb = 1 and (fib = 20°. In our computation of branching ratio and asymmetry we use LCSR 
models for heavy to light form factors both at non-zero and zero momentum transfer [41]. 
Our choices of input parameters are summarized in the Appendix. 

In Table 2, we display the experimentally (BaBar, CLEO and Belle) measured branching 
ratios and the weighted averaged values for the B + — > (fiK* + and B° — > (fiK* . The theoretical 
predictions in the SM for two different form factor models, the LCSR and BSW models are 
given in Ref.[40]. 



6.1 LR(RL) mixing contributions to B — > (f)K* 

In this section we will study the effect of LR(RL) mixing in B — > (fiK* process. This LR(RL) 
mass insertion can enhance the Wilson coefficient Cg g by a factor of nig /nib compared to the 
standard model contribution in the same way as shown in section 5 for B — > (fiK process. 
Hence, one need to impose the constrain on LR(RL) mixing from experimentally measured 
B(B — > X s 7) and also from S^k, C^k and B(B — > (fiK) as obtained section 5. 

In this scenario, the new weak phase ifi, (the phase of the LR(RL) mixing ) will contribute 
to direct CP-violating asymmetry Acp defined as : 

T{B+ ^(fiK*+)-T{B~ ^(fiK*-) 
CP T(B+ -> (fiK*+) + T(B~ -> (fiK*-) { } 
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X 


M) 


£ SUSY (in units of 10~ 6 ) 


A CP (in %) 


0.3 


(0.55 x 10- 2 , 1.8) 


31.83±i;g (32.45l|I) 


19.3918& (16.1718:8?) 


(0.82 x lO- 2 , 2.8) 


14.5011?! (21.6211^) 


20.2 (10.9618:8s) 


(0.82 x 10" 2 , 2.9) 


12.3911;^ (19.82±f;gf) 


15.73 (8.04l8:8i) 



Table 4: B(B+ -> <^*+) and „4 C p(5 + -> 0if*+) at the QCD scale p = m b /2 for LR mass 
insertion for selected points in the allowed p — ip space. The numbers in the parenthesis 
correspond to the RL mass insertion. The standard model branching ratio corresponding to 
this scale is (14.92l|?|) x 10~ 6 . The errors are due to ±10% theoretical uncertainties in the 
calculation. 

in terms of partial widths. Recently BaBar and Belle Collaboration has presented their 
measurement of CP violating asymmetries for B° — > <f)K*° and B^ — > (pK*^ 1 [43, 45] 

A C p(B° -> (f>K*°) = 0.04 ±0.12 ±0.02, 0.07 ± 0.1518:8a (39) 
A C p(B ± -> QK** 1 ) = ±0.16 ±0.17 ±0.04, - 0.13 ± 0.2918:?? (40) 

where, in each asymmetry result, the first number correspond to the BaBar data while the 
second one correspond to Belle measurement. The standard model value for this asymmetry 
is less than 1%. The new physics (SUSY) contributions from the new penguin operator 
appeared due to LR(RL) mixing (S b ^ RL ^) can modify the sign and magnitude of Acp(B ± — > 
(pK*^) within the allowed parameter space of S^ R ^ RL y 

To get the numerical values of B{B + -> <f)K* + ), and A C p{B + -> <f)K* + ), we fix x = 0.3 
and 4.0. Then for a given QCD scale //, we select some points in the allowed parameter space 
of p — if) plane (as marked by Z in Figure 2) for both values of x. In this computation, we 
include (±10%) theoretical uncertainties arising mainly from two sources. First, the different 
form factors, whose numerical values are highly model dependent. In our analysis we consider 
LCSR model [41], instead of BSW model [42]. Secondly, while computing the hadronic 
matrix element using the QCD improved factorization method, one encounters linear and 
logarithmic divergent terms. These divergent integrals are regularized by introducing some 
phenomenological parameters as explained earlier. However, these divergent integrals appear 
only in the transverse helicity amplitudes, whose contributions are small. Allowing for ±10% 
variation in the input values estimates the probable error due to these sources. 

In Table 3, we present the branching ratio B(B + — > <f)K* + ) and the CP rate asymmetry 
A(B + — > (j)K* + ) for p = rrib and selected values of x — 0.3 and 4.0 for values of p and tp 
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allowed by B — > data for LR mass insertion (RL is shown in the parenthesis). The 
branching ratio with SUSY turn out to be much higher than the standard model value of 
6.18li;i5, which is lower than the experimental data (Table 2). Even the lower range of 
theory prediction is much higher than the upper range of experimental data within la. The 
rate asymmetry has much less error and is consistently within the range ~ —4% to ~ 18%. 

Similarly in Table 4, we show B and Acp calculated for QCD scale p = m b /2. In this 
case, there are two allowed regions from the combined B — > <J)K and B — > X s ^y constraints 
corresponding to x = 0.3. For x = 4.0, there are no allowed regions from B — > <pK data. The 
standard model branching ratio is much larger compared to the one computed at p = mf,. 
In SUSY, apart from the QCD scale p, the branching ratio also depend on the values of p 
and ifi. Moreover, the selected points in p — ip plane are different in the two cases. In the 
case, with p = rrib/2, and at p = 0.82 x 10~ 2 and ip = 2.8,2.9 radians, with LR mixing, 
lower ranges of the theory predictions are consistent with the upper range of experimental 
data at one sigma. For the other value of p and ip, the theoretical prediction for branching 
ratio is much higher than the standard model theory as well as experimental data. With RL 
mixing, the predicted branching ratio is much larger compared to both the standard model 
prediction and experimental data. The asymmetries for both LR and RL mixing case are 
always positive with less errors. 

We conclude that for some selected points in p — ip plane allowed by B — > <pK and 
B — > X s 7 at p — m b /2 provide a satisfactory understanding of B — > <pK* process. We also 
note that, at p = mb the SUSY contribution to the branching ratio of B — > 4>K* is too large 
to be consistent with the experimental data. 

We have also studied other CP violating asymmetries that can arise in vector- vector final 
state. The set of observables are defined in terms of A Q , A» and A± as follows [30]. 



\A X \ 2 + \A X \ 2 „ \A X \*-\A > 



2 2 
A ±l = - Im (A ± A* - A ± A*) , A|| = Re + A\\A* ) , 

E ±i = -Im (A ± A* + A±A*^ , S|| = Re (A\\A^ - A\\A 



q 



P±l = Re p A* ± Ai + A*A ± , p ±± = lml ^-A* ± A 



Q 



q 



PHo = -Im (J [A^Ao + AlA^j , Pll = -Im ^Aj (41) 

where A = {0, ||, _L} and the observables where % = {0, ||}, We restrict ourselves to the 
study of helicity dependent CP asymmetry defined as Saa/Aaa [30]. For the purpose of 
illustration we select last two sample points from the Table 4. At these values of p and ip, 
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X 




^oo/ ^-00 


E||/A N 


S ±± /A ±± 


0.3 


(0.82 x 10~ 2 , 2.8) 


0.19 ±0.00 


o.6ilffi 


u -° ' -0.013 


(0.82 x 10- 2 , 2.9) 


0.15 ±0.00 


n 7 i +0.020 
u - ' x -0.026 


r, epr+0.022 
u - oo -0.027 



Table 5: Helicity dependent CP asymmetry at the QCD scale \x = rrib/2 for LR mass 
insertion for selected points in the allowed p—ip space. The errors consist of ±10% theoretical 
uncertainties. 

with LR mass insertion, the lower range of the B SUSY is consistent with the upper range 
of the experimental data on B[B — > 4>K*) at one sigma. We then compute Eaa/Aaa f° r 
each values of A for these two sets of p and ip and is shown in Table 5. As before, in this 
case also we include ±10% theoretical uncertainties in our calculation. We only show the 
helicity dependent asymmetries for LR mass insertion, since with RL mass insertion the 
SUSY contribution to the branching ratio is too large to be consistent with the data. 

7 Conclusions 

In this paper, we study the SUSY contribution to the gluonic dipole moment operator to 
B — > <pKs process. We find that the LR(RL) mass insertion can enhance the gluonic 
dipole moment operator significantly. We then use the experimentally measured quantities, 
such as S^k, C^k and B(B — > <pKs) to constrain the parameter space of LR(RL) mixing. 
Interestingly, we find that the constraints from B — > <pK data is consistent with the B — > X s 7 
limit. It turns out that the same enhancement of gluonic dipole moment operator can also 
affect other penguin dominated process, such as B — > <f)K*, which is a pure penguin process 
like B —> <fiK s . In standard model, the predicted Acp(B — > (j)K*) is less than 1%, which is 
consistent with the measured Acp- In SUSY, with LR(RL) mass insertion, one can have a 
new source of CP violation arising from the complex LR(RL) mixing parameter, which can 
provide additional contribution to the Acp{B — > <pK*). We calculate such asymmetries and 
also the branching ratio for the set of parameters allowed by B — > <pK data. At /i — m b , for 
both LR and RL mass insertion, we observe that the predicted branching ratio is well above 
the experimentally measured one. However, the Acp{B — > 4>K*) is consistent with the data. 
On the other hand, at the QCD scale \x = m b /2 with LR mass insertion, we find that the 
theoretically computed branching ratio is consistent with the data with in one sigma error. 
At this second choice of allowed parameter space of 8 b £ R , RL)) we find Acp(B + — > <pK* + ) 
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in the range 15% to 20%, which is significantly higher than the standard model prediction 
but is still consistent with the present data. Finally, we also present helicity dependent CP 
asymmetries in the same parameter space of S^ R . 
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8 APPENDIX : Input parameters and different form 
factors 

In this Appendix we list all the input parameters, decay constants and form factors used for 
the calculation of B — > <pK and B — > <pK* . 

1. Coupling constants and masses ( in units of GeV ): 

a em = 1/129, a s (m z ) = 0.118, G F = 1.16639 x 10~ 5 (GeV) -2 , 
a s (m z ) = 0.118, m z = 91.19, m b = 4.88, m B = 5.2787, 
m = 1.019, m K = 0.493, m K *± = 0.892 

2. Wolfenstein parameters : 

A = 0.2205, A = 0.815, r] = 0.324, p = 0.224, 

3. Constituent quark masses m^i = u, d, s, c, d)( in units of GeV): 
m u = 0.2, m d = 0.2, m s = 0.5, m c = 1.5, m b = 4.88. 

4. The decay constants (in units of GeV): 

/b = 0.19, U = 0.237, = 0.16, f K . = 0.221 

5. The form factors at zero momentum transfer : 

F^ K = 0.33, A l = 0.331, A 2 = 0.283, V = 0.458 
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The g 2 -dependence of A 1: A 2 and V in the LCSR models is given by 



A2{q 2 ) 



A^Oyi - 0.6(g 2 /m|) - 0.23(g 4 /™!) 
A 2 (0)/1 - 1.18(g 2 /m|) + 0.281(g 4 /™!) 
V(0)/1 - 1.55(g 2 /m|) + 0.575(g 4 /m|) 



(42) 
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Figure 1: Variation of | G(x) | with x(— m|/m?). 
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Figure 2: Contour plots of S^k, C^k and B(B — > <fiKs) in p — ip plane for two values of 
x = 0.3 (a, c) and 4.0 (b, d) for LR(RL) mixing with rriq = 500 GeV. The scale fi = mb for 
Figures (a) and (b), while it is m^/2 for Figures (c) and (d). The lcr allowed regions of S^k, 
B and C^k are two gray bands, area within the solid curves and area outside the two dotted 
contours respectively. 
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